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Abstract. 

Using new sets of stellar evolution models at very low metallicities (Z = 10 -7 , 10 -5 ) and previously published 
grids we examine spectral properties of the ionising continua, the Lyman-break, and the Lya and He n A1640 
recombination lines in starbursts. The metallicity dependence of these properties, especially the transition from 
primordial galaxies (Population III) to currently observed metallicities, is examined for various IMFs and star 
formation histories. 

For the average properties of starbursts, approximated by a model with constant star formation, the main findings 
are: 

— The Lyman continuum flux Q(H) increases with decreasing metallicity. For a universal Salpeter IMF from 1-100 Mq the 
enhancement reaches typically a factor of ~ 3 between solar metallicity and Pop III objects. 

— While for metallicities Z <; 1/50^© the amplitude of the Lyman-break depends little on Z, a reduction by a factor ~ 2 is 
found at lower metallicities, due to the strong increase of the average stellar temperature. 

— Using theoretical models and empirical constraints we discuss the expected evolution of the hardness of He + to H ionising 
photons, Q(He + )/Q(H), with metallicity and possible uncertainties. Over the metallicity range Z = to ~ 10 -4 the hardness 
decreases from log(Q(He + )/Q(H)) ~ — 1.4... — 2.3 by ~ 1.5-2 or more orders of magnitude, depending strongly on the upper 
mass cut-off of the IMF. From empirical constraints we derive a hardness log(Q(He + )/Q(H)) ~ -3.2 to -2.6 for metal-poor 
starbursts (1/25 ^ Z/Zq ^ 1/4) and softer spectra for higher metallicities. We also provide a simple estimate of the possible 
impact of hot WR like stars on Q(He + )/Q(H) at very low metallicities (Z ^ 10 -4 ). 

— Calibrations for star formation rate determinations from various recombination lines at all metallicities and for various IMFs 
are derived. 

For young bursts the maximum Lya equivalent width is shown to increase strongly with decreasing metallicity 
from VF(Lya) ~ 250-350 A at Z £ 1/50 Z Q to 400-850 A or higher at Z between 10" 5 and (Pop III) for the 
same Salpeter IMF. However, for well known reasons, the Lya emission predicted likely represents an upper limit. 
Non-negligible He II A1640 emission due to stellar photoionisation appears to be limited to very small metallicities 
(log(Z/Zo) —5.3) and Population III objects. 

should be 



The predictions, available on the Web through the CDS and at tittp : //webast . ast . obs-mip . f r/sf r/ 



useful for a variety of studies regarding high redshift galaxies, cosmological reionisation, and others. 

Key words. Cosmology: early Universe - Galaxies: stellar content - Stars: general - Stars: fundamental parameters 
- Stars: atmospheres 



1. Introduction 

The discovery of numerous high redshift galaxies provides 
a unique opportunity to study galaxies in formation in 
the early Universe. Most of these galaxies show signs of 
actively ongoing massive star formation, as revealed by 
their overall spectral appearance, by detailed spectral fea- 
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tuxes, and in many cases by the presence of strong emission 
lines. 

In fact, since the detection of Lyman-break galaxies at 
z ~ 2-4 by colour selection techniques (Steidel et al. 1996, 
review by Stern & Spinrad 1999), Lya surveys or other 
search techniques have found a large number of objects at 
higher redshift showing in most cases intense line emission 
(e.g. Hu et al. 1998, 1999, Kudritzki et al. 2000, Rhoads k 
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Malhotra 2001, Malhotra & Rhoads 2002, Ellis et al. 2002, 
Frye et al. 2002, Ajiki et al. 2002). This also includes the 
most distant galaxy known to date, a lensed galaxy at 
z = 6.56 found through its Lya emission (Hu et al. 2002). 

It is possible that such strong line emitters showing 
also relatively little continuum light may represent the 
earliest stages of galaxy formation, where small amounts 
of metals have so far been formed. Strong ongoing star 
formation and a small dust content, which can suppress 
Lya emission, would then explain the high observed Lya 
equivalent widths (cf. Hu et al. 1998). More striking is the 
suggestion of Malhotra & Rhoads (2002) that the high 
Lya equivalent widths observed in the LALA survey at 
z = 4.5 could, among other explanations, be due to metal- 
free (so-called Population III, hereafter Pop III) objects. 
Possibly we are beginning to probe distant chemically lit- 
tle evolved galaxies, closing the gap between the first (pri- 
mordial) galaxies and the high (close to solar) metallicities 
of massive galaxies in the local Universe. 

To properly study these objects appropriate spec- 
tral models are necessary. They should take into account 
all possible metallicities, and also probable systematic 
changes of the stellar initial mass function (IMF) at very 
low metallicity (Abel et al. 1998, Bromm et al. 1999, 
Nakamura & Umemura 2001, Hernandez & Ferrara 2001). 
Providing such model calculations is the main aim of the 
present work. 

Other applications also require an understanding of 
how properties like line emission and the ionising fluxes of 
starburst behave in the transition between metal-free (Pop 
III) objects and metal-poor galaxies with observable coun- 
terparts in the local Universe. This is the case in studies 
addressing the re-ionisation history of the Universe (e.g. 
Gncdin 1998, Ciardi et al. 2000, review by Loeb & Barkana 
2000), especially if account is taken for the simultaneous 
metal-enrichment (cf. Gnedin & Ostriker 1997, Ferrara & 
Schaerer 2002). Also, for searches of primordial galaxies, 
it is of interest to explore how far extreme properties such 
as strong He n emission predicted for Pop III starbursts 
(Tumlinson & Shull 2000, Tumlinson et al. 2001, 2002; Oh 
et al. 2001, Bromm et al. 2001b, Schaerer 2002) are truly 
limited to zero metallicity. Our model calculations, includ- 
ing in particular metallicities Z = 0, 10 -7 , 10 -5 , 4. 10 -4 
and higher, allow, for the first time, such investigations. 

The present paper is structured as follows. Our mod- 
els ingredients, including two new sets of stellar evolution 
tracks at very low metallicity, are described in Sect. ^. The 
predicted Lyman continuum fluxes and the properties of 
the Lyman-break at all metallicities are presented in Sects. 
H and|| respectively. In Sect. || we discuss theoretical pre- 
dictions and empirical constraints on the He + ionising flux 
and the hardness of the ionising spectra of starbursts at 
various metallicities. Finally, quantitative predictions for 
the Lya and He II A1640 emission are given in Sect. |. 
Section H summarises our main conclusions. 



2. Model ingredients 

The basic model ingredients are identical to those de- 
scribed in Schaerer & Vacca (1998, hereafter SV98) and 
the Pop III models of Schaerer (2002a, henceforth S02). 
A brief summary is provided subsequently including the 
new features introduced in the present work. 

2.1. Atmosphere models 

Depending on the metallicity different sets of atmosphere 
models are used. 

For metallicities Z < 10~ 5 we follow S02 in using 
the grid of plane parallel non-LTE atmospheres computed 
with the TLUSTY code of Hubeny & Lanz (1995) for 
T off > 20000 K and the plane parallel line blanketed LTE 
models of Kurucz (1991) with a very metal-poor compo- 
sition ([Fe/H] = — 5.) otherwise. Possible uncertainties in 
the predicted ionising spectra of very metal-poor stars are 
discussed in S02. In comparison with the computations of 
S02, the use of an extended grid of TLUSTY models in 
the present paper leads to some small differences related 
to the highest energy range considered here. 

For higher metallicities we follow SV98 and adopt for O 
stars the CoStar non-LTE models including stellar winds 
(Schaerer & de Koter 1997), for Wolf-Rayet (WR) stars 
the pure He spherically expanding non-LTE models of 
Schmutz et al. (1992), and Kurucz (1991) models of ap- 
propriate metallicity otherwise. 

As discussed in earlier publications (e.g. Mihalas 1978, 
Schaerer & de Koter 1997, Schmutz et al. 1992, Tumlinson 
& Shull 2000, Kudritzki 2002) the inclusion of non-LTE 
models is the most important ingredient to obtain accu- 
rate predictions for the ionising spectra of massive stars. 

2.2. Stellar tracks 

To explore a wide range of metallicities covering popula- 
tions from zero metallicity (Pop III) , over low metallicities 
(Z ~ 4. x 10~ 4 ) such as observed in H n regions in the local 
Universe, up to solar metallicity (Z = Zq = 0.02), we use 
the following stellar evolution tracks: 1 ) the Pop III tracks 
covering masses from 1 to 500 M Q with no/negligible mass 
loss compiled in S02 (Marigo et al. 2001, Feijoo 1999, 
Desjacques 2000), 2) new main sequence stellar evolution 
tracks from 1 to 500 M© computed with the Geneva stel- 
lar evolution code for Z = 10 -7 , 3) non-rotating stellar 
models from 2 to 60 M Q from Meynet & Maeder (2002) 
complemented with new calculations for 1 M Q and 85-500 
Mq for Z = 10~ 5 , 4) non-rotating Geneva stellar evolu- 
tion tracks for masses 0.8-120 M Q (up to 150 M Q for 
Z = 0.0004) from the compilation of Lejeune & Schaerer 
(2001) for the remaining metallicities Z =0.0004, 0.001, 
0.004, 0.008, 0.02 (=Z ), and 0.04. For massive stars we 
have adopted the high mass loss tracks in all cases, as these 
models reproduce best various observational constraints 
from the local Universe (cf. Maeder & Meynet 1994). 
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Table 1. Line emission coefficients q in [erg] for Case B, 
n e = 100 cm" 3 and the different adopted T e . See S98 and 
S02 for sources of the atomic data 



Line 



Cl 

T e =30 kK 



T e =10 kK 



appropriate Qi 



Lya 

He ii A1640 
Ha 



1.04 x 10" 11 

5.67 x 10" 12 
1.21 x 10" 12 



1.04 x 10" 11 

6.40 x 10~ 12 
1.37 x 10" 12 



Q(H) 

g(He+) 

Q(H) 



Our calculations at Z < 10~ 5 include only the H- 
burning phase. As He-burning is typically less than 10 
% of the main sequence lifetime, and is generally spent 
at cooler temperature, neglecting this phase should have 
little or no consequences on our predictions. A possible 
exception may be if stars at very low Z are very rapid 
rotators (cf. Maeder et al. 1999) which could suffer from 
non-negligible rotationally enhanced mass loss and could 
therefore become hot WR-like stars (cf. Sect. |6.2| ). 

To verify our new computations we have compared sev- 
eral Z = 10~ 5 tracks with the independent calculations 
of Meynet & Maeder (2002) done with a strongly mod- 
ified version of the Geneva stellar evolution code. Good 
agreement is found regarding the zero age main sequences 
(ZAMS), H-burning lifetimes, and the overall appearance 
of the tracks. 

The HR-diagram showing the main sequence tracks at 
very low metallicity is given in Fig. [I]. As expected one 
finds an important shift of the ZAMS and main sequence 
toward hotter T c g with decreasing Z. For massive stars 
the trend shown by the low Z tracks is expected to con- 
tinue down to a limiting metallicity Zu m of the order of 
10~ 9 , below which the stellar properties essentially con- 
verge to those of metal-free (Pop III) objects. This is the 
case, as massive stars (M ^ 50 M Q ) with < Z < Zn m 
will rapidly build up a mass fraction Acno ~ 5. 10 -10 to 
2. 10~ 9 of CNO during pre main sequence contraction or 
the early main sequence phase (e.g. El Eid et al. 1983, 
Marigo et al. 2001). A typical value of Z\i m = 10 -9 will be 
adopted subsequently for various simplified models fits. 



2.3. Evolutionary synthesis models 

As described in S02 the above stellar atmosphere models 
and evolutionary tracks have been included in the evolu- 
tionary synthesis code of Schaerer & Vacca (1998). Using 
the prescriptions summarised below we compute the pre- 
dicted properties of integrated stellar populations at dif- 
ferent metallicities for instantaneous bursts and constant 
star formation, the two limiting cases of star formation 
histories. 




Fig. 1. HR-diagram showing the main sequence tracks of 
stars with masses ^ 3 to 500 (150 for Z = 0.0004) M Q 
of metallicities Z = (Pop III, dotted line), Z = 10~ 7 
(solid), Z = IQ- 5 (long-dashed), and Z = 4. 10~ 4 (short- 
dashed) and their ZAMS. The source of the tracks is given 
in Sect. ^. Note the important shift of the ZAMS to high 
T ff from low metallicity to Z = 0. 



2.3.1. Nebular emission 

Among other predictions the code in particular computes 
the recombination line spectrum including Lya, He n 
A1640, He n A3203, He i A4026, He i A4471, He n A4686, 
H/3, He I A5016, He I A5876, and Ha. Case B recombina- 
tion is assumed for an electron temperature of T e = 30000 
K at Z < 10~ 5 and T e = 10000 K otherwise, and a low 
electron density (n e = 100 cm -3 ). Lya emission is com- 
puted assuming a fraction of 0.68 of photons converted 
to Lya (Spitzer 1978). The line emission coefficients q 
(defined by Eq. |t]) of interest here are listed in Table [l]. 

Some differences, e.g. due to a more realistic temper- 
ature structure or due to collisional effects on H lines, 
can be expected between the adopted prescriptions and 
predictions from detailed photoionisation models (see e.g. 
Stasihska & Tylenda 1986, Stasihska & Schaerer 1999). 
However, for the scope of the present investigation these 
effects can quite safely be neglected. 

As shown by S02 the inclusion of nebular continuous 
emission processes is crucial for very metal-poor objects 
with intense ongoing star formation. Following S02 we in- 
clude free-free, free-bound, and H two-photon continuum 
emission assuming T e = 20000 K or 10000 K and the above 
value of n e . 
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Fig. 2. Predicted SEDs including Lya and He n A1640 emission lines for zero age main sequence (ZAMS) models at 
different metallicities. Left panel: Dependence of the SED on metallicity: The metallicities Z = 0. (Pop III), 1CP 7 , 
10~ 5 , 0.0004, and 0.02 (solar) are from top to bottom in the EUV (A < 912 A), and reversed at longer wavelengths. 
The dashed lines are the pure stellar emission, the solid lines show the total (stellar + nebular) emission. IMF A in all 
cases. Right panel: Dependence of the SED on M up for metallicities Z = 0. (Pop III), 10 -7 , and 10~ 5 . Solid lines show 
the IMF A, dashed lines the IMF B. In all cases the total (stellar + nebular) emission is shown. Discussion in text 



Table 2. Summary of IMF model parameters. All models 
assume a Salpeter slope for the IMF 



Model ID 


M low 


M up 


cm 


A 


1 


100 


2.55 


B 


1 


500 


2.30 


C 


50 


500 


14.5 



2.3.2. Initial mass function 

In view of our ignorance on the massive star IMF at very 
low metallicities (Z < 4. x 10 -4 ) we adopt as in S02 a 
powerlaw IMF and different upper and lower mass limits 
with the aim of assessing their impact on the properties 
of integrated stellar populations. The main cases modeled 
here are summarised in Table || []. For all models the IMF 
slope is taken as the Salpeter value (a = 2.35) between 
the lower and upper mass cut-off values Mi ow and M up 
respectively. 

The model A IMF is a good description of the IMF in 
observed starbursts (e.g. Leitherer 1998, Schaerer 2002b) 

1 For comparison with SFR indicators used for "normal" 
galaxies we also indicate the mass conversion factor cm ex- 
pressing the relative masses between our model IMF and a 
"standard" Salpeter IMF with Mi ow = 0.1 and M up = 100 
Mq used frequently throughout the literature, cm is given by 
cm = J" M$(M)dm/ J^ up M$(M)dm. 



down to - 1/50 Z Q (= 4. 10~ 4 ), the metallicity of I Zw 18 
representing the most metal-poor galaxy known to date. 
It is adopted in all calculations for Z > 4. 10~ 4 . IMFs 
B and C, favouring the formation of very massive stars, 
could be representative of stellar populations at metallic- 
ities Z ^ Z cx it ~ 10 -5 (e.g. Bromm et al. 2001a), where 
altered fragmentation properties may form preferentially 
more massive stars (cf. Abel et al. 1998, Bromm et al. 
1999, Nakamura & Umemura 2001). The computations at 
Z < 10~ 5 consider all IMF cases (A, B, and C). Note that 
our calculations obviously do not apply to cases where 
only one or few massive stars form within a pre-galactic 
halo, as suggested by the simulations of Abel et al. (2002). 

3. Overview of EUV-UV spectra of starbursts at 
various metallicities 

To illustrate several points discussed in detail below, we 
plot the EUV to UV spectral energy distribution (SED) 
for selected models at different metallicities and computed 
with different IMFs (see Fig. |). Only ZAMS models are 
shown here for simplicity sake. 

The main points apparent from this Figure are: 

— The rapid decrease of the He + ionising photon flux 
(i.e. flux above 228 A) and the associated overall soft- 
ening of the ionising spectrum from Pop III to higher 
metallicity (cf. Sects. [| and ||). 
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Fig. 3. Temporal evolution of the H ionising photon 
flux Q(H) (upper panel) and the hardness <2(He + )/<2(H) 
(lower panel) for instantaneous bursts at all metallicitics Z 
between 0.04 (=2 Z Q ) and 0. (Pop III). The very metal- 
poor models (Z = 10~ 5 , 10~ 7 and 0) with the IMF B 
(1-500 M Q ) are shown as solid lines from top to bot- 
tom. The remaining metallicities computed for the IMF A 
(1-100 M ) are shown with dashed lines. Their hardness 
(2(He + )/Q(H) reached for constant star formation at equi- 
librium is shown by the solid triangles (Z < 10 -5 ), and 
the short lines for larger Z. Predictions for IMF cases not 
shown here are summarised in Tables || and ^. Discussion 
in text 

— The increase of the Lyman-break amplitude with in- 
creasing Z (Sect. ||). 

— The rapid decrease of the contribution of nebular con- 
tinuous emission to the UV continuum flux (cf. S02). 

— The strong dependence of the ionising flux and in turn 
also of the nebular continuous emission on the upper 
mass cut-off of the IMF, which is considered as a free 
model parameter for very low (Z < 10 -5 , currently 
unobserved) metallicities (Sects. |]aiid^). 

The metallicity dependence of the emission lines 
(strongest in the Pop III models) cannot clearly be dis- 
cerned from this figure. It is discussed in Sect. |^. 

4. Ionising properties of starbursts at various 
metallicities 

4.1. Burst models 

The basic quantities describing the ionising spectrum are 
the emitted number of H, He, and He + ionising photons, 
denoted by <2(H), <2(He°) and <2(He + ) respectively, and 
the hardness Q(He+)/Q(H) (Q(He°)/Q(H)) tracing the 



Table 3. ZAMS properties of integrated stellar popula- 
tions at various metallicities. All quantities are normalised 
to a total burst mass of 1 M . The logarithm of the ion- 
ising photon fluxes Qi in units of [photon s _1 M _1 ] is 
given. 



7 

Zj 


TMP 

11V1.T 


Q(H) Q(He°) 
[log( (photon s~ 


^/(M© yr 




n 
u. 


A 

1\ 


46.98 


46.75 


40.04 


A P. 90 
40. ZZ 


n 
u. 


R 


47.29 


47.10 


AP OP 
40. ZO 


AP. /in 
40. 4U 


n 
u. 


p 


47.98 


47.80 


ai n^ 

4 f .UO 


AP QP 

40. yo 


in- 7 

1U 


\ 


46.94 


46.65 


A Q A K 
4O.40 


AP IP 
40. OO 


in- 7 


D 
13 


A 7 Qfl 

4/.oU 


A 7 f\« 
4/.UD 


A ^ £1 
40.01 


AP KK 
40.00 


1U 


P 

\J 


48.01 


47.78 


AP QQ 

40. oy 


A 1 1 A 
4 i . 14 


10~ 5 


A 


46.90 


46.55 


42.39 


46.44 


1(T 5 


B 


47.30 


46.99 


44.56 


46.64 


1(T 5 


C 


48.02 


47.73 


45.35 


47.24 


0.0004 


A 


46.88 


46.43 


43.4 


46.45 


0.0004 


A* 


47.01 


46.58 


43.62 


46.53 


0.001 


A 


46.86 


46.41 


43.44 


46.47 


0.004 


A 


46.85 


46.38 


43.41 


46.50 


0.008 


A 


46.84 


46.35 


43.67 


46.52 


0.020 


A 


46.85 


46.34 


43.72 


46.54 


0.040 


A 


46.87 


46.33 


43.71 


46.59 



*: Salpeter IMF with M up =150 M and Mi ow =l M Q 

energy range between 54 (24.6) and 13.6 eV. The pre- 
dicted temporal evolution of <2(H) is shown in Fig. ^ (up- 
per panel) for all metallicities between Z = 0. and 2 x Zq. 
For the very low metallicities (Z < 10 -5 ) only the models 
with an IMF extending to 500 M (model B) are shown 
for clarity sake. Adopting a larger value of M up affects 
only the predictions at very young ages (ages ^2.5 Myr) 
due to very short lifetime of these stars. 

The predicted Qi of ZAMS populations (age=0) for 
all IMF cases are listed in Table |[ For completeness with 
S02 the photon flux in the Lyman- Werner band (11.2-13.6 
eV) capable to dissociate H2 is also listed. 

The main difference in the Lyman continuum photon 
output at different Z is a slower decline of the ionising 
photon production at low metallicities, due to the blue- 
ward shift of the main sequence. The temporal evolution 
of Q(H) at Z = 10~ 7 is essentially undistinguishable from 
the Pop III case. The larger <2(H) apparent for Z < 10 -5 
at ages ^2.5 Myr are essentially due to the larger value of 
M up adopted at very low Z. The difference at older ages 
(when stars with masses > 100 M have disappeared) 
represents the pure metallicity difference. 

As can be seen from Table || the <2(H) production 
of ZAMS populations at different metallicities increases 
somewhat with decreasing Z; the changes remain fairly 
small (Js 40 %) in reasonable agreement with other esti- 
mates (e.g. Tumlinson & Shull 2000). However, in cases 
such as constant star formation (equivalent to a temporal 
average) the Z-dependence is more pronounced (cf. Sect. 

O). 
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4.2. Constant star formation 

The main predictions for models with constant star for- 
mation at all metallicities and for all the IMF cases are 
listed in Table |^. In this case most quantities of interest 
here reach rapidly (over timescales Si 6-10 Myr; except 
for the Lyman-break and Q(LW) requiring i> 200 Myr) 
an equilibrium value given in the table, normalised to a 
star formation rate (SFR) of 1 M Q yr -1 . In addition to 
the ionising photon production Qi (cols. 3-5), and the Hi 
photodissociating photon flux (Q(LW), col. 6), we list the 
average energies E(Q(R)) and E(Q(Hc + )) of the Lyman 
continuum photons and the photons with energies above 
54 eV (cols. 7 and 8). These quantities, not further dis- 
cussed here, are e.g. of interest to estimate the thermal 
evolution of the ISM. Most of the data for Z = 0, 0.0004, 
an 0.02 were already given in Table 3 of S02. Due to the 
use of a finer grid of atmosphere models at Z = some 
small changes are found for these models^. The values in 
Table ||] supersede those of S02. 

As expected from the earlier discussion (see Fig. |J), 
the Lyman continuum flux Q(H) shows an increase with 
decreasing metallicity, which can be fitted to an accuracy 
better than 10 % by 

log(Q(H)) = -0.0029 x (log(Z) + 9.) 2 ' 5 + 53.81 (1) 

valid for the Salpeter IMF from 1-100 M Q (models A) and 
Z G [10~ 9 , 0.04]. The limitation to Z > Z Vlm ~ 10~ 9 is due 
to the fact that all metallicities below Z\\ m are equivalent 
to the Pop III (Z = 0) discussed in Sect. ^2[ Fits 

for alternate IMFs may be derived from the data in Table 

i 

Overall, while ZAMS Lyman continuum fluxes vary by 
less than <, 40 % over the entire metallicity range for IMF 
A (Table ||) , the ionising output at SFR= const shows an 
increase of an factor ~ 1.9 (2.8) between solar and 1/50 
Zq (zero metallicity). Even larger increases are of course 
predicted in the case of IMFs extending to higher masses 
(models B and C). 

5. Lyman-break predictions at various metallicities 

Quantitative predictions of the amplitude of the Lyman- 
break, i.e. the relative flux above and below the Lyman 
edge are of importance for example to estimate the es- 
cape fraction / csc of ionising photons out of their hosts. 
Such measurements are possible in relatively nearby [z 
0.015) and high redshift galaxies (e.g. Leitherer et al. 1995, 
Deharveng et al. 1997, 2001, Steidel et al. 2001). 

Different working definitions of the Lyman-break ex- 
ist in the literature. For simplicity we chose the definition 
adopted in Starburst99 (Leitherer et al. 1999), quantifying 
the break amplitude by 912+/912 - , where 912 + is the av- 
erage flux in F\ units over the interval 1080-1120 A, and 
912~ the average over 870-900 A. A word of caution is 

2 /i640 and the corresponding Q(He + ) are reduced by ~ 9, 
4, and 30 % for IMFs A, B, and C respectively. 
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Fig. 4. Temporal evolution of the Lyman-break for instan- 
taneous bursts and constant star formation at all metal- 
licities < Z < Zq . Dashed lines show the predictions for 
Z > 4. 10~ 4 , solid lines are for Z = (Pop III), Z = 10~ 7 , 
and 10 -5 (from bottom to top). The behaviour of the very 
low Z burst models at ages logt ~ 7.2, 7.5, 7.7 are artifacts 
due to the limited numerical resolution in tracks. Note the 
small sensitivity of the Lyman-break for SFR= const at 
metallicities Z > 4. 10~ 4 and the important decrease at 
lower metallicities. 

appropriate if wavelengths A ^ A (Lya) are used as refer- 
ence for the red side of the Lyman-break. In this case the 
continuum may, for very metal-poor objects, include a non 
negligible contribution from nebular continuous emission 
(see S02) which itself also depends on the escape fraction 
of Lyman continuum photons. A detailed SED fit should 
then be undertaken to quantify / esc . 

The temporal evolution of the Lyman-break predicted 
by the models at metallicities < Z < Z & is shown in Fig. 
^ for instantaneous bursts and constant star formation. 
The main finding apparent here is the overall reduction 
of the Lyman-break at metallicities Z 10 -4 " , e.g. by 
0.2-0.3 dex at very young ages, for identical IMFs or larger 
values of M up . This is due to the decrease of the break 
in individual stars with increasing T c ff and to the shift of 
the main sequence towards hotter temperatures. At higher 
metallicity (Z ^ 4. 10 -4 ) this temperature/metallicity de- 
pendence is not important, and e.g. the predicted values 
at SFR= const (near equilibrium at t i> 10 s yr) vary by 
less than 10 % around log(912+/912 - ) ~ 0.58. In con- 
trast for very low metallicities one has for the same case 
log(912+/912 _ ) ~ 0.3-0.4 for the IMFs A and B. 

3 IMFs A and B yield Lyman-break amplitudes differing by 
less than ~ 0.1 dex. 
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Table 4. Predictions for the case of constant star formation models (equilibrium values) with different IMFs (col. 2) 
at all metallicities (col. 1). Given are the ionising photon fluxes Qi (cols. 3-6), their average energies E (cols. 7, 8), 
and the proportionality constants fi between recombination line luminosities and the SFR (cf. Eq. 0, cols. 9-11). The 
logarithm of the ionising photon fluxes Qi is given in units of [photon s _1 M Q -1 ]. 
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*: Salpeter IMF with M up =150 Mq and Mi ow =l M 

uncertain predictions for reasons discussed in Sect. ^ 
": no data provided, as predictions uncertain and strongly overestimated (cf. Sect. 



Note that, as expected, at Z > 10~ 3 our predictions 
at young ages are somewhat (~ 0.1 dex) smaller than 
the Starburst99 models for identical Z and IMF, due 
to the inclusion of non-LTE O star model atmospheres 
in the present computations. For constant star forma- 
tion this difference becomes, however, smaller. We also 
compared our calculations to the models of Smith et al. 
(2002; Norris, private communication) including more so- 
phisticated non-LTE stellar atmospheres. For metallicities 
0.001 < Z < Zq and SFR= const their Lyman-break pre- 
dictions show an additional but small (~ 10-20 %) reduc- 
tion. 



6. The Q(Re + )/Q(R) hardness of starbursts at 
various metallicities 

6.1. Predictions at metallicities Z = to Z Q 

6.1.1. Burst models 

The time evolution of the predicted hardness 
Q(He + )/Q(H) is shown in the lower panel of Fig. 0. 
The corresponding hardness reached at equilibrium in 
the case of constant star formation (cf. below), plotted 
on the right, illustrates the non negligible difference with 
the hardness derived from a simple ZAMS population 
neglecting stellar evolution effects (cf. S02). 

Note that predicted quantities such as Q(Hc + )/Q(H), 
rely obviously strongly on the adopted value of M up . For 



the very low Z this quantitative dependence can be esti- 
mated from the tabulated ZAMS properties (Table ||). 

As expected from the strong decrease of the stel- 
lar temperatures with increasing metallicity (cf. Fig. ^) 
both the maximum hardness (at age=0) and the overall 
Q(He+)/Q(H) decreases for Z between and 10~ 5 . The 
typical timescale for a decrease of Q(Hc + )/Q(H) by ~ 2 
dex in a burst is driven by the redward stellar evolution, 
and is short (~ 2-3 Myr), with obvious potential implica- 
tions for the detection of sources with very hard spectra 
(cf. Sect. |?|). Possible additional sources of He + ion ising 



photons not included here are discussed in Sect. p. 2 . 

At higher metallicities (Z > 4. x 10 -4 ) the present 
models predict a re-increase of Q(He + )/Q(H) after J> 3-4 
Myr, due to presence of WR stars, among which a frac- 
tion is found at high temperatures (cf. Schmutz et al. 1992, 
SV98). Albeit with minor quantitative differences, a qual- 
itatively similar behaviour is predicted by the Starburst99 
models based on very similar input physics (Leitherer et 
al. 1999). However, these predictions depend especially on 
the procedure adopted to link stellar tracks with atmo- 
spheres in WR phases with strong winds, and on the ne- 
glect of line blanketing in the adopted WR model atmo- 
spheres. The reality and the extent of such a trend remains 
therefore questionable, especially at the largest metallici- 
ties (cf. review of Schaerer 2000). Indeed using recent line 
blanketed O and WR atmospheres and different prescrip- 
tions to connect the interior and atmosphere models Smith 
et al. (2002) find a considerably softer spectrum - i.e. re- 
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Fig. 5. Hardness Q(He + )/Q(H) of the He + ionising flux 
for constant star formation as a function of metallicity (in 
mass fraction) for all models given in Table || At metal- 
licities above Z > 4. 10~ 4 the predictions from our models 
(crosses), as well as those of Leitherer et al. (1999, open 
circles), and Smith et al. (2002, open triangles) are plot- 
ted. The shaded area and the upper limit (at higher Z) 
indicates the range of the empirical hardness estimated 
from H II region observations (see Sect. 3.3). At very low 
Z solid lines show the fits to the data given by Eq. 0. 
Discussion in text 

duced Q(He + )/Q(H) - before and during the WR phase. 
To circumvent this theoretical uncertainty we will subse- 
quently derive an empirical estimate of the hardness at 
metallicities Z k, 4. x 10" 4 (Sect. O). 



6.1.2. Constant star formation 

For the case of constant star formation at equilibrium 
the metallicity dependence of the hardness Q(He+)/Q(H) 
of the ionising flux is shown in Fig. || for all the IMFs 
considered. As apparent, for the very metal-poor cases 
(Z & 10~ 4 ) Q(He+)/Q(H) can be well fitted by 



log(Q(Hc+)/Q(H)) = a x log(Z) + b 



(2) 



with the numerical coefficients listed in Table |5|, again tak- 
ing into account that all metallicities Z < Zn m ~ 10 -9 are 
equivalent. 

At higher metallicities — in the Z range of known ob- 
jects — the theoretical predictions for Q (Hc + ) / Q (H) are 
probably less clear, due to possible presence of hot WR 
stars, difficulties in their modeling (cf. above), and the 
neglect of non-stellar emission processes which could con- 
tribute especially to Q(He + ). Indeed as shown in Fig. ^| 
rather important differences are obtained between vari- 



Fig. 6. Additional contribution Q(He + )*/Q(H) of puta- 
tive hot WR-like stars to the total He + hardness of the 
ionising flux Q(He+)/Q(H) (cf. Eq. g) for SFR= const 
as a function of the minimum mass Mwr. Shown are the 
cases for iwR.A* = 10 and 50 % (solid and dashed curves), 
temperatures T = 100 and 160 kK (lower and upper curve 
respectively), and M up =500 (100) M Q extending over the 
full plot (lower left only). 



ous evolutionary synthesis models (SV98, Starburst99 of 
Leitherer et al. 1999, and the latest computations of Smith 
et al. 2002). In the metallicity range Z ~ 8. 10~ 4 to 5. 10~ 3 
our empirical estimate of Q(He + )/Q(H) (Sect. |6.3| ), also 
shown on this Fig., is likely more reliable than the mod- 
els. At still larger Z the average ionising spectrum of star- 
bursts should be softer, as indicated by the tentative em- 
pirical upper limit and predicted by the Smith et al. (2002) 
models. 



6.2. Uncertainties on the He + ionising flux at very low 
metallicities (Z <, 10~ A ) 

In contrast to the Lyman continuum flux Q(H), the He + 
ionising flux (and in general spectral features at high en- 
ergy) show a very strong dependence on the stellar tem- 
perature in the T e g range ~ 50-100 kK typical of very low 
metallicity stars (Fig. 2 in S02). Therefore their prediction 
is naturally more sensitive to even small modifications of 
the exact stellar T e ff or evolutionary scenario. 

For example, one may wonder how reliable the above 
predictions of the metallicity dependence of Q (Hc + ) /Q (H) 
(Fig. ||) are at Z 10~ 4 , where presently no observational 
constraints are available. In fact, studies of massive stars 
in the Local Group suggest that their average rotation 
rates increase towards low Z (Maeder et al. 1999), which 
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- when combined with their increased compactness - can 
lead to non-negligible mass loss despite the low metallic- 
ity (Maeder & Meynet 2000, Meynet & Maeder 2002). If 
this effect is large enough, one could imagine that fast 
rotators could loose sufficient mass to follow a WR star 
like evolution leading possibly to a He/C/O core at tem- 
peratures T <; 100 kK, a scenario known for metal-rich 
massive WR (e.g. Maeder & Meynet 1988). Despite high 
rotational velocities, the detailed calculations of Meynet 
& Maeder (2002) for Z = 10~ 5 do not show important 
alterations of the evolution for stars with M < 60 M©. 
Exploratory calculations of Marigo et al. (2002) for Pop 
III stars treating in an simplified manner the effects of 
rotation on stellar mass loss find such a scenario for stars 
with initial mass ^ 750 M Q . 

Quantitatively the effect of such a putative hot "WR- 
like" population on the hardness of the ionising flux can be 
estimated for the case of constant star formation only in 
the following way. Suppose that stars of given initial mass 
M > Mwr spend this phase at constant luminosity L and 
(hot) temperature T during a constant fraction /wr = 
(^wrA*) of their lifetime t*. Assuming their winds are 
optically thin at > 54 eV, the He + flux in this phase is 
then 



Q(He+)*(M) = ^M 9 (He+r(T), 



(3) 



where q(Re + )* is the photon flux per unit surface area 
(cf. Fig. 2 in S02). The hardness (Q(He+)/Q(H))' of a 
population including both the "normal" (main sequence) 
population and the additional hot population can be ap- 
proximated by (for SFR= const) 



/ Q(Hc+) V Q(He+ 
^ Q(H) ) ~ Q(H) 



(1 - /wr) + 



Q(He+) 
Q(H) 



WR, 



(4) 



assuming small changes in the total Lyman continuum 
photon output Q(H). Here the first term stands for the 
unperturbed "normal" hardness shown in Fig. ^. The sec- 
ond term describes the additional contribution due to hot 
WR-like objects, which can be expressed as 



Q(He+)* _ £ U(M)Q(Ke+y(M)$(M)dM 
Q(H) 



U{M)Q{K){M)${M)dM 



(5) 



where Q(H)(M) is the lifetime average of the Lyman con- 
tinuum production, and $(M) is the IMF. The quantities 
t*(M) and Q(H)(M) are taken from S02, and we assume 
as first approximation a luminosity L(M) corresponding 
to the ZAMSQ. The hardness contribution Q(He+)*/Q(H) 
is then computed assuming temperatures T = 100-160 kK 
as expected from stellar evolution models, and durations 
/wr = (*WrA*) °f 0-1 or 0-5- The result is plotted for 
the IMFs A and B (M up =100 and 500 M ) in Fig. | as 
a function of the minimum mass A/wr, above which all 
stars are assumed to reach this hot "WR" phase. 

4 The results discussed below are fairly insensitive to the 
exact (low) metallicity adopted for these quantities. 



Table 5. Fit coefficients for Eq. || 
IMF a b 



A 
B 

C 



-0.66 ± 0.071 
-0.37 ± 0.028 
-0.39 ± 0.028 



-8.22 ± 0.51 
-5.04 ± 0.20 
-4.98 ± 0.20 



The case of (£wrA*) — 50 %, which would require 
very strong mass loss already during the main sequence or 
a nearly homogeneous evolution leading early to a blue- 
ward evolution, appears extremely unlikely and is shown 
here to mimic the "strong mass loss" models adopted in 
the Pop III models of S02. For a hot phase of a duration 
typical of the post main sequence evolution (~ 10 % of 
total lifetime) Fig. [| shows that such putative "hot WR" 
could in the "best" case contribute an additional hardness 
Q(He + )*/Q(H) of the order of 10~ 2 - -3 , comparable to 
the hardness of normal stellar populations with metallici- 
ties Z ^ 10~ 4 (cf. Fig. |5|). To examine how realistic such 
cases may be, will require a detailed understanding of the 
coupled processes of stellar mass loss, rotation, and inter- 
nal mixing. At present the available limits are Mwr > 60 
M at Z = 10~ 5 and (*wr/**) ^ 10 % from the rotating 
stellar models of Meynet & Maeder (2002), and M WR > 
750 Mq at Z = from the simplified models of Marigo et 
al. (2002). 

Although the above exercise shows that at very low 
metallicity (Z < 10~ 5 ) the hardness Q(He+)/Q(H) due 
to stellar sources could be higher than shown in Fig. 0, it 
seems that such scenarios are quite unlikely. If star forma- 
tion takes place on much longer time scales, and massive 
stars would not form (or in much smaller quantities), hot 
planetary nebulae could also be a source of hard ionis- 
ing photons, as illustrated by the scenario of Shioya et al. 
(2002). In any case, a major uncertainty stems from our 
limited knowledge of the IMF at very low metallicities. 

6.3. Empirical constraints on the He + ionising flux of 
starbursts at Z ^ 4. x 10" 4 

Spectroscopic observations of extra-galactic giant H n re- 
gions probing He II recombinations lines can yield empiri- 
cal information on the "average" hardness Q(He + )/Q(H) 
of starbursts. Indeed it is well known that a fairly large 
fraction of metal-poor H n regions show the presence of 
nebular He n A4686 emission indicative of a hard ionis- 
ing spectrum (see e.g. Guseva et al. 2000, compilation of 
Schaerer et al. 1999). A complete explanation of the ori- 
gin of the required high energy photons (shocks, X-rays, 
WR stars) remains to be found (e.g. Garnett et al. 1991; 
Schaerer 1996, 1998; Guseva et al. 2000, Izotov et al., 2001, 
Stasinska & Izotov 2002). 

The largest sample of high quality data is that of Izotov 
and collaborators (cf. Guseva et al. 2000 and references 
therein), which shows He n A4686 detections with typical 
relative intensities of J(4686)/7(H/3) ~ 1-2 %. 
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From such a sample we may estimate an average hard- 
ness from 



Q(He+) _ 7(4686) 
Q(H) ~ J(H/3) 



x a x /detect x 



tde 



ttotal 



(6) 



where o (=0.47 (0.63) for T e = 10000 (30000) K) trans- 
lates the relative line emissivities, /detect is the fraction 
of objects showing nebular He n among the spectra of 
sufficient quality to allow detections down to this inten- 
sity level, tdetect is the age range covered by the observed 
H II regions, and ftotai is the total lifetime of H n re- 
gions. Inspection of the data of Izotov yields /detect ~ 1/3- 
1/2 for metallicities Z G [8. 10~ 4 , 4.8 10" 3 ] and possibly 
even a larger fraction at the lowest observed metallicities. 
(Stasihska & Izotov 2002). One has ttotal ~ 10 Myr, the 
typical H n region lifetime or the time over which bursts 
show detectable line emission. However, it is now well es- 
tablished that H II region samples suffer from selection bi- 
ases leading to an absence of advanced bursts, or in other 
words a preferential selection of bursts with ages younger 
than < 4-5 Myr (e.g. Bresolin et al. 1999, Raimann et 
al. 2000, Moy et al. 2001, Stasihska et al. 2001). For this 
range of values and a = 0.47 one obtains an estimate of 
log(Q(He+)/Q(H)) 3.2 to -2.6. Postulating the ex- 
treme case that additionally the first 2 Myr of the H II 
phase are also not detected in the optical^ one has a min- 
imum value of ^detect ~ 1 Myr, yielding a lower limit of 
log(<5(He + )/Q(H)) <; -3.5. In principle this average value 
could be even larger as the (conservative) assumption of 
^detectAtotai < 1 implies that older (non-detected) H II 
regions have no He + ionising flux. 

This estimate is obviously independent of the nature 
of the hard (He + ) ionising radiation. By construction Eq. 
^ provides an estimate of the average Q(He + )/Q(H) ex- 
pected in objects with constant ongoing star formation for 
metallicities Z ~ 8. 10~ 4 to 5. 10 -3 . The near absence of 
nebular He n detections in H II regions at higher metallic- 
ity (cf. Schaerer 1998, Guseva et al. 2000) indicates softer 
spectra. However, it is difficult to establish a firm upper 
limit on Q(He+)/Q(H) for Z £ 5. 10~ 3 . We here retain 
log(Q(He + )/Q(H)) < -3.3 as a tentative limit. 

6.4. Overall behaviour of Q(He+)/Q(H) with 
metallicity 

In short, from the considerations above, we find the follow- 
ing two cases for the most plausible metallicity dependence 
of the average Q(He + )/Q(H) hardness ratio of starbursts 
with metallicity (see Fig. ||). 

1) If a universal Salpeter like IMF with a "normal" 
upper mass limit of ~ 100 M prevails for all metallicities 
the hardness decreases by more than 2 orders of magnitude 
from Z = to ~ 10 -4 , re-increases thereafter (up to Z ^ 
5. 10~ 3 , in metal-poor starbursts) to a level ~ 2 to 10 



5 The well known absence of H II regions with very large H/3 
equivalent widths could be taken as an indications for such a 
scenario. 



times smaller than that of Pop III objects, and decreases 
again to low levels for higher metallicities. 

2) If very massive stars are favoured at metallicities 
Z < 10~ 4 , the hardness Q(He+)/Q(H) of Pop III objects 
is considerably enhanced (corresponding to a powerlaw 
spectrum with spectral index a ~ 2.3-2.8 in F v ), then 
decreases down to levels somewhat smaller or comparable 
to that of metal-poor starbursts, before decreasing further 
to levels at least two orders of magnitude softer than at 
zero metallicity. 

7. Lya and He n A1640 emission at various 
metallicities 

To first order recombination line luminosities are propor- 
tional to the ionising photon flux and are thus simply ex- 
pressed as (in units of erg s^ 1 ) 



Li{t) [ergs" 1 ] = ci (1 - / csc ) Q t {t) [s~ 



(7) 



where q are the line emission coefficients given in Table 
0, / csc is the photon escape fraction out of the galaxy or 
observed aperture, and Qi is the ionising photon flux (in 
units of s -1 ) corresponding to the appropriate recombina- 
tion line. 

However, as well known Lya constitutes a particular 
case due to its very large line optical depth, which implies 
that several effects (e.g. dust absorption, ISM geometry 
and velocity structure) can alter the total Lya emission 
and lead to complex line profiles (cf. Chariot & Fall 1993, 
Valls-Gabaud 1993, Chen & Neufeld 1994, Kunth et al. 
1998, Loeb & Rybicki 1999, Tenorio-Tagle et al. 1999). 
Furthermore for Lya source at redshifts close to or above 
the redshift of reionisation the intrinsic Lya emission may 
be further reduced or suppressed by absorption in nearby 
line of sight HI clouds (cf. Miralda-Escude & Rees 1997, 
but also Haiman 2002, Madau 2002). A proper treatment 
of these effects requires a complex solution of radiation 
transfer which depends strongly on geometrical properties 
of the ISM and IGM, and for which no general solution 
is possible. One must thus caution that depending on the 
application our simplifying assumptions may not apply 
and the predicted Lya emission should thus be considered 
as an upper limit. Note that these difficulties do not affect 
other recombination lines such as Ha and He II A1640, 
whose optical depth is strongly reduced compared to Lya. 

Bearing the above in mind, the time evolution of the 
Lya and He n A1640 line luminosities can be deduced from 
the evolution of Q(H) and Q(He + ) respectively given in 
Fig. § 

In the case of constant star formation, at equilibrium, 
recombination line luminosities L/ are proportional to the 
star formation rate (SFR), i.e. 



Ll = (1 - /csc)// 



/ SFR 



(8) 



For Li in erg s 1 the proportionality constants fi for Lya, 
Ha, and He n A1640 are listed in columns 9-11 of Table 
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Fig. 7. Temporal evolution of the Lya equivalent width (left panel) and He II A1640 equivalent width (right) for 
instantaneous bursts at all metallicities. The very metal-poor models (Z = 0, 1(T 7 , and 1(T 5 ) with the IMFs C 
(50-500 M Q ), B (1-500 M Q ) and A (1-100 M Q ) are shown as short-dashed, solid, and dotted lines respectively from 
top to bottom. The remaining metallicities (for IMF A) are shown with dashed lines. The equilibrium values for 
SFR= const at metallicities Z < 10 -5 are plotted on the right (at arbitrary ages) using open squares for the IMF C, 
filled triangles for IMF B, open circles for IMF A, and using short lines for higher metallicities (with IMF A). Note 
the very large maximum W(Lya) predicted at young ages. IF(He n A1640)^ 5 A are only expected at the lowest 
metallicities [Z ^ 10~ 7 ), e xcept if hot WR-like stars not included in the tracks were formed e.g. through important 
stellar mass loss (Sect. 3.2). See text for further discussion 



||. The variations of the ionising flux Q(H) already dis- 
cussed above imply in particular lower star formation rates 
at low metallicity when identical recombination lines are 
used. We note that our Ha SFR conversion factors are in 
good agreement with other computations at various metal- 
licities (e.g. Kennicutt 1998, Sullivan et al. 2001) when 
rescaled for their Mi ow using cm (Table ||). 

The predicted Lya and He II A1640 emission line equiv- 
alent widths of ageing bursts of different metallicities and 
all the IMF cases considered are shown in Fig. 0. Note 
that our revised Pop III models show smaller Lya equiv- 
alent widths compared to the previous calculation in S02. 
This is due to an erroneous continuum definition in the 
earlier computations. The new models supersede those of 
S02. Good agreement is also found with the IF (Lya) pre- 
dictions of Tumlinson et al. (2002). The reader is also 
reminded that none of these recent calculations include 
stellar Lya absorption (cf. Valls-Gabaud 1993, Chariot & 
Fall 1993, and Valls-Gabaud & Schaerer 2002 for new pre- 
dictions) . 

Maximum Lya equivalent widths of ~ 240-350 A are 
predicted for metallicities between solar and 4. 10~ 4 . For 
Z down to zero (Pop III), max(IF(Lya)) may reach up 
to - 800-1500 A for the various adopted IMFs (cf. S02) 
For comparison, the equilibrium values for SFR=const are 
in the range of IF(Lya) - 175-275, 240-350, 500-930 A 



for IMF A, B, and C respectively at Z < 10~ 5 , and ~ 
70-100 A for higher metallicities (IMF A). Note that the 
increased Lyman continuum output of young very metal- 
poor populations alone does not explain the strong in- 
crease of IF (Lya) (cf. Table |). In addition the reduced 
stellar UV continuous luminosity at A ~ 1215 A, due to 
the shift of the SED peak far into the Lyman continuum 
(Fig. ^), also contributes to increase W(Lya). 

A high median Lya equivalent width (~ 430 A) was 
found in the Large Area Lyman Alpha (LALA) survey 
of Malhotra & Rhoads (2002) at z = 4.5 and interpreted 
as due to AGN, starbursts with flat IMFs, or even Pop 
III objects. Indeed, if constant star formation is appro- 
priate for their objects and the IMF slope is universally 
that of Salpeter, the observed large IF(Lya) would require 
very metal-poor populations with large upper mass cut- 
offs and/or an increased lower cut-off (e.g. IMFs B or C). 
Alternatively their observations could also be explained 
by predominantly young bursts, with metallicities Ss 10 -5 
and no need for extreme IMFs (Fig. This issue will 
be addressed in detail in a subsequent publication (Valls- 
Gabaud & Schaerer 2002). 

As expected from the softening of the radiation field 
with increasing metallicity, the He n A 1640 equiva- 
lent widths decreases strongly with Z; values TF(He n 
A1640) <; 5 A are expected only at metallicities Z ^ 10~ 7 , 
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except if hot WR-like stars (cf. Sect. 3.2) or non-stellar 
sources (e.g. X-rays, AGN) provide sufficient amounts of 
He + ionising photons. 

Part of the relative weakness of H^He n A1640) com- 
pared to W(Lya) is due to a relatively strong, mostly 
nebular, continuum flux at 1640 A (see S02). As Lya emis- 
sion may be strongly reduced due to the effects discussed 
earlier in objects beyond the re-ionisation redshift, and 
the He II A1640 luminosity is potentially strong enough 
to be detected (cf. Tumlinson et al. 2001, Oh et al. 2001, 
Schaerer & Pello 2001), it is a priori not clear if both lines 
may be observed simultaneously and if so which of the two 
lines would be stronger. 



8. Conclusions 

We have examined the spectral properties of the ionising 
continua, the Lyman-break, and the Lya and He II A1640 
recombination lines in starbursts with metallicities Z from 
zero - corresponding to primordial, Pop III objects - over 
low metallicities [Z <; 4. 10~ 4 ) observed in nearby galax- 
ies, up to solar metallicity (Z =O.O2). 

Our computations, including new sets of stellar evolu- 
tion models at very low metallicities (Z = 10~ 7 , 10 -5 ) and 
previously published grids at other allow us in particu- 
lar to study how spectral properties vary in the transition 
from Pop III objects to "normal" currently measurable 
metallicities. 

Various IMFs are treated, including also cases where 
very massive stars (up to ~ 500 M Q ) are formed, as sug- 
gested by hydrodynamical calculations for metallicities 
Z ^ Z clit w 10~ 5 (Bromm et al. 2001a, Abel et al. 1998, 
Nakamura & Umemura 2001). 

Predictions are provided for the number of H, He , 
and He + ionising photons and average photon energies in 
these continua, the hardness of the ionising spectrum, the 
amplitude of the Lyman-break, the number of photons 
able to photodissociate H 2 , and finally recombination line 
luminosities and equivalent widths (mostly for Lya and 
He ii A1640). 

Two limiting cases of star formation histories, instan- 
taneous bursts and constant star formation, are consid- 
ered. For SFR=consi, presumably appropriate to describe 
the average properties of starbursts galaxies or popula- 
tions thereof, the following main results have been ob- 
tained: 

• As expected from numerous earlier computations the 
Lyman continuum flux Q(H) increases with decreasing 
metallicity. For a universal Salpeter IMF from 1-100 
Mq the enhancement reaches typically a factor of ~ 
3 between solar metallicity and Pop III objects for a 
constant star formation (Sect. ^, cf. Tumlinson & Shull 
2000, Schaerer 2002). 

• While for metallicities Z J; 1/50Z Q the amplitude of 
the Lyman-break is rather metallicity independent, a 
reduction by a factor ~ 2 is found at lower metallicities, 



related to the strong increase of the average stellar 
temperature (Sect. ||). 

• The predicted hardness of the ionising fluxes between 
> 54 eV and > 13.6 eV, i.e. the ability to doubly ionise 
He, decreases by ~ 1.5-2 or more orders of magnitude 
from Z = to ~ 10~ 4 depending strongly on the upper 
mass cut-off of the IMF (Sect. |). 

From empirical constraints we derive a hardness 
logQ(He + )/Q(H) ~ -3.2 to -2.6 for metal-poor star- 
bursts (1/25 ^ Z/Zq ^ 1/4) and softer spectra for 
higher metallicities (Sect. |6.3| ). The former should pro- 
vide the best estimate of Q(He + )/Q(H); the latter 
finding is also compatible with recent evolutionary syn- 
thesis models of Smith et al. (2002) including line blan- 
keted non-LTE atmosp here s for WR and O stars. 
We also provide (Sect. |6.2[ ) a simple estimate of the 
possible impact of hot WR like stars on Q(Hc+)/<5(H) 
at very low metallicities (Z 10~ 4 ). Such stars could 
eventually form due to a strong enhancement of mass 
loss related to rapid rotation (Marigo et al. 2002) or 
in principle also due very efficient rotational mixing 
processes (cf. Meynet & Maeder 2002), although we 
believe that these scenarios are quite unlikely or overall 
of minor importance. 

• Finally, calibrations for star formation rate determi- 
nations from Lya and other recombination lines at all 
metallicities and for various IMFs are provided (Sect. 

0). 

For young bursts, the maximum Lya equivalent width 
predicted is shown to increase strongly with decreasing 
metallicity from W(Lya) ~ 240-350 A at Z ^ 1/50 Z Q to 
400-850 A or higher at Z between 10~ 5 and (Pop III) 
for the same Salpeter IMF (Sect. @). 

We find that non-negligible He II A1640 emission due 
to photoionisation from stellar sources appears to be lim- 
ited to very small metallicities (log(Z/Z Q ) ^ —5.3) and 
Population III objects, except if hot WR like stars, whose 
existence appears very speculative, were frequent. 

The detailed model predictions presented here are 
available on the Web through the CDS and at 
tittp : //webast . ast . obs-mip . f r/sf r/. In subsequent 



publications our models are applied to the interpreta- 
tion of Lya observations in high redshift galaxies (Valls- 
Gabaud & Schaerer 2002), modeling of the combined 
chemical enrichment and re-ionisation history of the 
Universe (Ferrara & Schaerer 2002), and feasibility stud- 
ies on the detection of Population III objects (Pello & 
Schaerer 2002). 
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